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We investigate the QCD phase diagram of isospin asymmetric matter using the Polyakov loop extended quark 
meson (PQM) model with vector interaction. The critical point temperature is found to decrease in isospin 
asymmetric matter and disappear at large isospin chemical potential. We also discuss the QCD phase transition 
in the neutron star core. From comparison of the QCD phase diagram in PQM and corresponding baryon and 
isospin chemical potentials of neutron star matter in relativistic mean field models, we show that the order of the 
chiral phase transition in the neutron star core could be crossover because of large isospin chemical potential. 
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I. INTRODUCTION 

The QCD phase transition would be realized not only in 
heavy-ion collisions but also in compact astrophysical ob- 
jects and phenomena such as heavy neutron stars |ll|l, super- 
novae [21] and black hole formations Jj-la]. At zero baryon 
chemical potential (/ie), the QCD phase transition at finite 
temperature (T) is accessible by using the lattice Monte Carlo 
simulation, e.g. 13, S]. At large fig , fis/T > 1, the situa- 
tion is much less clear since the lattice simulation is plagued 
by the well-known sign problem 191]. We can investigate this 
region by using chiral effective models such as the Nambu- 
Jona-Lasinio (NJL) model lUOll and the quark meson (QM) 
model ][lly, and the chiral effective model with the Polyakov 
loop effects such as the Polyakov loop extended Nambu- 
Jona-Lasinio (PNJL) model Ill2l - [l4ll and the Polyakov loop 
extended quark meson (PQM) model iH [H]. The QCD 
phase diagram, especially the QCD critical point (CP) loca- 
tion, strongly depends on models and model parameters IUTIi . 
Therefore, further experimental and theoretical developments 
are necessary to determine the structure of the QCD phase di- 
agram. 

For laboratory experiments, the search for CP in heavy-ion 
collisions is ongoing at RHIC I ISO and is planned in the com- 
ing FAIR facility. Since the phase transition is second order 
at CP, the coherence length ^ is divergent and large fluctua- 
tions of the order parameter are expected in a volume of the 
size <^'^. Various signatures of CP have been proposed the- 
oretically ill9ll . It is not an easy task to observe the diver- 
gence signature of ^ in heavy-ion collisions, since the system 
size and the evolution time are limited. Moreover, it is diffi- 
cult to create cold dense matter, and CP may not be reachable 
in the laboratory if CP is located in the high density region, 
Atfl > 500 MeV. 

By comparison, very dense matter is formed in compact 
astrophysical phenomena. For example, high density and 
low temperature matter is formed in the neutron star core, 
and high temperature and high density matter is produced 
during a gravitational collapse of a massive star and binary 
stars 1201 . From the observation of these phenomena, we 



may get information on the QCD phase diagram in the high 
density region (ll-lg]. In compact astrophysical phenomena, 
charge neutrality leads to suppressed proton fraction com- 
pared with that of neutrons, and the isospin chemical potential 
Sfi = {fjLn ~ f^p)/'^ = if^d — fJ-u)/'2 is finite and positive. In 
particular, Si^t appears as another independent thermodynam- 
ical variable in supemovae and BH formations, since trapped 
neutrinos modify the neutrinoless charge neutrality condition 
{Sfi = Me/2). Therefore it is necessary to consider (5/i depen- 
dence of the QCD phase diagram in order to discuss the QCD 
phase transition in compact star phenomena. 

The phase structure in the three thermodynamic variables 
(T, /i, (5/i) is still an open problem. In our previous work la, 
we have discussed the possibility of the CP sweep during BH 
formation processes where Sfi is finite; quark matter core and 
hadronic envelope may merge to one phase, when the tem- 
perature exceeds the CP temperature (Tcp). The location 
of CP strongly depends on Spi; for large 6fi, Tcp becomes 
lower and it becomes more probable for the heated matter 
to go through CP. There are several recent works which dis- 
cuss the QCD phase diagram in charge neutral dense mat- 
ter 1I21I1 and in the three-dimensional space, (T, /i, S^) ll22ll23ll 
or (T, /^, /Xi) 1I24I1 . where fi^ is the lepton-number chemical 
potential. The phase diagram structures in these works have 
some differences. In Ref. 12 ill, the isospin chemical poten- 
tial is found to be small Sfi < m7r/2 in charge neutral quark 
matter, and pions are not found to condense. In Ref. ]122|]. 
three-dimensional (T, fi, Sfi) phase diagram is investigated in 
the mean-field approximation, and the s-wave pion condensed 
phase is found to appear in the finite /i and Sfi region. Tqp 
decreases with increasing Sf-i, until thejCP hits the pion con- 
densation phase boundary. In Ref. 112 3ll . fluctuation effects 
are taken into account in the quark meson model by using the 
functional renormalization group (FRG) flow equation. Tcp 
is also found to decrease with increasing Sfi. The s-wave pion 
condensed phase is found in the high Sfi and low fi region, 
but it is suppressed at large fi. As a result, the pion condensed 
phase is separated from the chiral first order phase transition 
surface in the (T, /x, Sfi) space. In Ref. ll24ll . Tcp is found to 
be insensitive to the lepton-number chemical potential. 



In this article, we investigate the isospin chemical potential 
dependence of the QCD phase diagram in more detail, and 
discuss the order of the chiral phase transition in the neutron 
star core, where T = and /iB,(5/i > 0. For this purpose, 
we first compute the QCD phase diagram using the two-flavor 
PQM with vector interaction and examine the 5^ dependence 
of the QCD phase diagram. According to the s-wave -kN re- 
pulsion argument [|25] and functional renormalization group 
results II23I1 , we assume that pions do not condensate. We then 
discuss the order of the chiral phase transition in neutron star 
core, through the comparison of the QCD phase diagram with 
the /3 equilibrium line in neutron star matter calculated by us- 
ing hadronic equation of states (EOSs). 

The article is organized as follows: In Sec. 11, we briefly 
describe PQM with vector interaction. The results are dis- 
cussed in Sec. in, where we show the dji dependence of the 
QCD phase diagram and compare the QCD phase diagram in 
PQM with neutron star matter chemical potentials. Sec. IV is 
devoted to summary and discussion. 



and P are the Polyakov loop and its conjugate. 



P = l^Tri , P = -r^Trit 



L is defined in the Euclidean space as. 



(5) 



L = P cxp i / dTA4 



(6) 



where P stands for the path ordering. The logarithmic term 
hi H{P, P) in Eq. (O comes from the Haar measure of the 
group integral in strong-coupling lattice QCD lUSll . Coeffi- 
cients, a{T) and h{T), are given as functions of T, and pa- 
rameterized as a{T) = ao + ai{Tn/T) + a2{To/T)^, and 
b{T) = h{T„/Tf Q. 



B. Effective Potential 



II. POLYAKOV LOOP EXTENDED QUARK MESON 
MODEL 

A. PQM Lagranglan and parameters 



We now give the effective potential in dense asymmetric 
matter in PQM. In asymmetric matter, u and d quark popula- 
tions are unbalanced, and we need to introduce two indepen- 
dent chemical potential for u and d quarks 



In this Section, we describe the PQM model augmented 
with the vector interaction. PQM is an effective model 
which has the chiral symmetry and confinement property of 
QCD in [ilL The Lagrangian density of the two-flavor PQM 
is given bylll[ll] 



C^q [ij^Df, - g{a + ij^T ■ tt) 



l^i9,,a 



gui^^ii 



hd,7rf-U{a,^) 



^U{P,P,T), 



gpl^r ■ H^] q 



R^ ■ Rn 



(1) 



fJ-u 



fi- Sfj. , fj.d = fJ- + Sfi , 



(7) 



where /i = /^te/S is the quark chemical potential. The isospin 
chemical potential Sfj, is an independent thermodynamical 
variable in supemovae or black hole formation processes, 
while the neutrino-less /3-equilibrium condition, 6fi = /ie/2, 
applies to cold neutron star matter 

We assume that the a meson and the temporal components 
of w and p'^ mesons take finite expectation values, while others 
do not. These expected values are assumed to be constant. In 
this approximation, the quark single-quasiparticle energy is 
given by 



where q denotes a quark field with Dirac, color and flavor in- 
dices, T is the Pauli matrix in the flavor space and o;^"^ and 
Rf^'^ are the field tensors of ui and p mesons. The mesonic 
potential U and the Polyakov loop potential U are given as. 



f/(cr, 7r) = A((t2 + tt^ - ^2)2/4 - ha 
.4/ «(r) 



U[P,P,T] =T' 



-PP + b{T)lnH{P,P) 



H{P, P) = l- 6PP + 4(P'^ + p3) - 3{PPf 



(2) 
(3) 
(4) 



a and tt are the isoscalar-scalar and isovector-pseudoscalar 
meson fields. The covariant derivative D^ = d^ — iA^ in 
Eq. ([T]i is the Dirac operator with a temporal static and homo- 
geneous background gluon field A^j = S^qAq. Without the 
explicit symmetry breaking term, the last term in Eq. (|2), the 
Lagrangian in Eq. ^ has SU{2)l x SU{2)ii symmetry. 
U{P, P, T) is an effective potential of gluon field, where P 



Z7* 



Ep + g^uj + gpT^R , 



with 



Ep = ^/p^ + AP, M ^ ga 



(8) 



(9) 



u! and R in Eq. (|8]l denote the expectation values of ui and p" 
mesons {ui ~ {ojo),R ~ (Rq))^ respectively, where the sub- 
script denotes the temporal component and the superscript 
for R shows isospin. The effect of vector interaction is to shift 
the quark chemical potential ll26ll . For later convenience, we 
define effective chemical potentials for u and d quarks, 

Pu = p-5p- gujuj - gpR , pd = P + Sp- gu,oj + gpR . 

(10) 

Integrating over the quark fields results in the following effec- 



tive potential, 

^PQM = U{P, P, T) + U{a, TT = 0) + Oo + f^T , (11) 
f7o = -2NjN,, Jj^EpO (A2 - p2) ^ (12) 



f]. 



1 



m^ cj2 + m2i?2) 






i^l = 1 + 3Per 

Fi^l + 3Pe-^< + 3Pe-2/5< + e-^^^'l , (15) 



£±^ Ep± jlf 



g-CS/3t_ ^ (14) 

(15) 
(16) 



where VIt is the thermal contribution and SIq is the fermion 
vacuum energy, regularized by the ultraviolet cutoff A. This 
term is necessary to reproduce the second-order chiral phase 
transition at zero baryon chemical potential jIb in the chiral 
limit |;16il . Each term on the right hand side of Eq. (fT4l i corre- 
sponds to the thermal contribution of zero, one, two, and three 
quark states. Similarly, Eq. (fT5t is the thermal contribution 
of antiquarks. While PQM is renormalizable and we can use 
dimensional renormalization lll6ll . it is sufficient to cut large 
momenta by a hard cutoff for our purposes. 

The equations of motion are obtained from the stationary 
conditions in equilibrium. 



d^ 


dn 


dn 


on 


dn 


da 


dp 


dP 


dbJ 


dR 



0. 



(17) 



We obtain (T, ^^ , 5^i) dependence of the mean fields, 
tr, P, P, (jj and R, by solving these equations. 

Here we do not consider the pion condensation because the 
s-wave pion condensation will not occur in neutron stars when 
we take account of the s-wave TriV repulsion Ii25i1 . Functional 
renormalization group analysis also shows the shrinkage of 
the pion condensed region at finite jjl than naively expected 
{5^JL > mJT) m. 



C. Model parametrization 

The parameter in the scalar-pseudoscalarpart, g, A, v^ h are 
fixed to reproduce some properties of quarks and mesons in 
vacuum for a given value of the hard momentum cut off A = 
600 MeV in this work. The quark-scalar meson coupling g 
is determined by the constituent quark mass in the vacuum 
niq = ga = 335 McV. The mesonic potential parameters A 
and V are given by the chiral condensate in the vacuum a = 
/tt = 92.4 MeV, and the a meson mass to^ = d^n/da^ = 
(700 MeV)^. The explicit symmetry breaking parameter h is 
given by the pion mass h = rn^fTj. 

In this study, we assume the quark-vector couplings are the 
same {g^ = Qp = 9v) for simplicity. We regard g,j as a 
free parameter, and we compare the results with several val- 
ues of r = gv/g- We also assume the common vector meson 
masses (?Ti^ = nip = my = 770 McV) 



The parameters in the Polyakov loop potential are fitted to 
the pure gauge lattice data luTll . The standard choice of the 
parameters reads |il4| oq = 3.51, ai = —2.47, 02 = 15.2 
and 63 = —1.75. The parameter Tq in Eq. (O sets the de- 
confinement scale in the pure gauge theory, i.e. Tq = 270 
MeV. Chemical potential dependence of these parameters is 
not considered in this work iflSl [2811 . 



III. RESULTS 

In this section, we discuss the Sfj. and the vector coupling 
dependence of the QCD phase diagram. The chiral phase tran- 
sition is found to be weakened at finite Sjj, or with finite vector 
coupling r. In order to demonstrate this point, we first discuss 
the order parameters as functions of yUg at several values of 
Sfi and r. 

The phase structure is obtained from the behavior of the or- 
der parameters a, P and P. Figure [T] shows /i^ dependence 
of the order parameters, a (left) and P (right), for several 
isospin chemical potentials at T = 96.5 MeV — Tcp{S^i = 
50 McV, r = 0) (CP temperature at S^ = 50 MeV and the 
vector-scalar coupling ratio r = 0). For small J/i, the chi- 
ral phase transition is first-order, while for S^ > 50MeV, the 
chiral phase transition becomes crossover The change of the 
nature of the phase transition with the increase of 6fi is not a 
peculiarity of the PQM model; in fact, several chiral models 
share this property, as discussed in Q (see also the Appendix 
for a discussion within the NJL model). 

In Fig. 121 we show a (left) and P (right) as functions of the 
baryon chemical potential at T = 101.5 MeV = Tcp((S^ = 
0, r = 0.2) and Sfj. = MeV for several values of the vector- 
scalar coupling ratio r. For the strong vector interaction, the 
transition chemical potential is shifted to higher values, and 
the chiral phase transition is smoothed. The transition be- 
comes crossover for r > 0.2 at this T. 

We next discuss the S/j. and vector coupling dependence of 
the chiral and deconfinement phase boundaries. Since the chi- 
ral phase transition at small i^i^ is actually a smooth crossover 
for finite quark masses, we have to establish a criterion to 
identify the phase boundary of the chiral transition. Here we 
define the chiral critical temperature Tc or the baryon chemi- 
cal potential ^b.c of the chiral phase transition by the peak of 
the chiral susceptibility Xa as a function of T or /ie for fixed 
Sfi and yUs or T, respectively. Since Xa- is divergent at the crit- 
ical point, we can unambiguously determine the critical point 
temperature Tqp and baryon chemical potential /icp by the 
diverging peak of Xa- in the T ~ jjlb plane. The chiral sus- 
ceptibility is defined as the second derivative of the effective 
potential by the explicit chiral breaking coefficient h. 



Xc 



= -T 



3d '{n/T) 



(18) 



Since h oc m^ is proportional to the bare quark mass, 
the above definition gives a susceptibility which is propor- 
tional to the usual definition around the critical point, Xa = 
—d'^n/dM^/T'^, where M is the bare quark mass. We nor- 
malize Eqs. ( fTSl ) by multiplying some powers of T to con- 
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FIG. 1. The order parameters cr(left panel), P(right panel) as functions of baryon chemical potential ^_b at T = 96.5 MeV = Tcp{Sn — 
50 MeV, r — 0) and three different isospin chemical potentials 5^ = (solid line), 50 (dash line), 70 MeV (dash-dot line). The vector-scalar 
coupling ratio is chosen to be r = 0. 
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FIG. 2. The order parameters a(left panel), P(right panel) as functions of baryon chemical potential /is at T = 101.5 MeV = Tcp((5/i 
MeV, r = 0.2) and 5/i = MeV for several values of the vector-scalar coupling ratio r = 0, 0.2, 0.3. 



sider dimensionless susceptibility. Figure |3] shows the chiral 
susceptibiUty as a function of temperature for several baryon 
chemical potentials at (5/i = McV. For each fig, we find 
a peak in Xtr, where the chiral phase transition occurs. At 
CP, {T,^b) = (TcpjMCp), this quantity is divergent which 
signals a second order phase transition. The critical points are 
found to be (Tcp, a^cp) = (117, 975) McV at (5/i = without 
vector coupling r = 0. 

As for the case of the chiral transition, the deconfinement 
transition is a crossover for finite quark masses, and we need 
to specify a criterion to identify the deconfinement phase 
boundary. Several prescriptions to define the critical tem- 
perature for deconfinement have been used in the literature: 
the temperature at which the Polyakov loop susceptibility, 
Xp, is maximum; the temperature at which dP/dT is max- 
imum OOl I31I1 ; finally, the half- value prescription, in which 
one identifies the deconfinement temperature with the average 
of the temperatures at which at P = 1/2 and P = 1/2 ||29I1 
(the two differ at finite /i). The Polyakov loop susceptibil- 
ity Xp may have a double peak structure in some cases lUSll : 
one peak is related to the chiral phase transition and the other 
is related to the transition caused by the Polyakov loop mean 



field potential. A similar double peak behavior is found in 
dP/dT 1I30II31I1 . Thus it is not easy to unambiguously define 
the deconfinement temperature from the Polyakov loop sus- 
ceptibility or the temperature derivative. Since the Polyakov 
loop is small (P, P ~ 0) in confined phase and large (P, P ~ 
1) in deconfined phase, the half-value prescription is the sim- 
plest one to adopt. Therefore, we adopt the half-value pre- 
scriptions to define the deconfinement temperature. 

We show the QCD phase boundaries for several 5^ values 
in Fig. m The hadron phase shrinks a little and the critical 
point temperature Pcp decreases with increasing 5^, while 
the confinement-deconfinement phase boundary only weakly 
depends on Sji. The reduction of the transition chemical po- 
tential may be understood as the density effects. For a simple 
estimate, let us consider the low T transition in the chiral limit 
without the vector coupling, where the sum of u and d quark 
number densities in the chiral restored phase is proportional 
to (yu + S^if + ifi- Sn)^ = 2^3(1 + iS^'^/n'^) as in the 
free massless case . If the QCD phase transition at finite (5/i 
occurs at the same density in the Wigner phase as that for 
(5/i = 0, the transition quark chemical potential is calculated 
to be /i ~ A'c — <5/i^//ic, where /ic represents the transition 
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FIG. 3. The chiral susceptibility Xct as a function of temperature for 
different isospin chemical potential /is — O(solid line), /icp(dash 
line), 1000 MeV(dash-dot line) at 5^ = MeV. x<t 's divergent at 
the critical end point. 



250 

200 

>150 

HlOO 

50 





5^=0,50,70MeV 



r=0 




200 400 600 800 1000 1200 

liB (MeV) 



FIG. 4. The QCD phase diagrams for several isospin chemical po- 
tential. The red dash and the blue solid lines show the crossover 
and first order chiral phase transition boundaries respectively at 
Sfj. — 0,50,70 MeV. The square dots show the CP. The black dash- 
dotted lines show the confinement-deconfinement phase boundaries 
at 6/1 = 0, 70MeV. 



chemical potential at J/i = 0. This estimate gives the transi- 
tion chemical potential shifts of 7.2 and 14 MeV for dfi = 50 
and 70 MeV, respectively, which is comparable to the PQM 
results, 7.0 and 13 MeV. Another possible explanation is the 
decrease of the effective number of flavors. At finite 5/i, one 
of the M or d quarks is favored, and the phase diagram is ex- 
pected to be closer to that at Nf = 1, where the phase transi- 
tion is weaker 

We note that the deconfinement transition temperature (Td) 
is a little higher than the chiral transition temperature Tc- The 
present behavior is consistent with the lattice Monte-Carlo 
simulation results, which suggest Td > Tc 13311 . It should 
be noted that this order depends on the choice of Tq and is 
different from some of the effective model results [32]. While 
the order of Td and Tc at /i = is an interesting problem on 
the relation of deconfinement and chiral transitions, it is irrel- 



evant to our conclusion and we choose Tq = 270 McV in the 
later discussion. 

It should be noted that the deconfinement phase boundary is 
almost insensitive to the baryon chemical potential, leading to 
a splitting of the chiral and deconfinement transition bound- 
aries. This behavior is similar to the strong coupling lattice 
QCD results including finite coupling and Polyakov loop ef- 
fects II31I1 . but it is different from the results obtained from 
the functional renormalization group method starting from the 
PQM initial condition at large cutoff Q. 
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FIG. 5. The QCD phase diagrams for several quark- vector meson 
couplings. The red dash and blue solid lines show the crossover 
and first order chiral phase transition boundaries respectively at 
r = 0, 0.2, 0.3. The black dash-dotted lines show the confinement- 
deconfinement phase boundaries at r = and r = 0.2. 

Figure|5]shows the QCD phase diagrams of symmetric mat- 
ter {6fj. ~ 0) for several quark-vector meson couplings. With 
increasing vector coupling, the chiral phase boundary moves 
to the higher fig direction, and the CP moves to the higher 
jig and lower T direction. The behavior of /icp is understood 
from the effective fis shift. We can ignore the p° meson ef- 
fects in symmetric matter, and the effective chemical potential 
is given as /2 = /i — rgoj. Therefore, a strong vector interac- 
tion makes fl small for a given ^^ BSll . and the phase bound- 
aries and the CP moves to high /ig for finite vector coupling, 
r ^ 0. By comparison, the vector coupling dependence of the 
confinement-deconfinement phase boundary is small. 

We show the QCD phase diagram in (T, /x^ , S/i) space in 
Fig. |6] As already mentioned, (5/i reduces Tqp and the tran- 
sition baryon chemical potential at T = 0. Then the first or- 
der boundary narrows with increasing 6fj,, and eventually the 
CP disappears at a certain value of 5/x. This happens also 
for the NJL model, as we discuss in more detail in the Ap- 
pendix. This behavior is important when we consider the chi- 
ral phase transition in dense and isospin asymmetric matter, 
which is realized in the core of neutron stars, where S/^i be- 
comes large. For example, the reduced CP temperature may 
affect the dynamical black hole (BH) formation processes. 
The highest temperature during the BH formation is calcu- 
lated to be T ^ 70 McV, and compressed matter may experi- 
ence either the first order, crossover, or CP sweep depending 
on the CP location in asymmetric matter [|6|]. 
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FIG. 6. The 3-dimensional (T, fis , S/i) QCD phase diagram in PQM 
with r — 0.2. The dash lines and the blue solid lines show the 
crossover and first order chiral phase boundary respectively. The red 
solid line shows CP for different S/i. 
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FIG. 7. This figure shows the first order phase boundaries at T = 
in PQM model with r = 0(red solid line), 0.2(red dash line) and 
0.3(red dot line) and {fis, Sfi) for neutron star matter calculated by 
NJL Q with r = 0(black dash-dot fine), TMl (Hcblue dot line) 
and lOSTY (Hcblue dash-dot fine) (H. 



Another example of dense asymmetric matter is the neu- 
tron star core. In the neutron star core, the internal tem- 
perature is of the order of 10® K r^ 10^^ MeV, which is 
small enough compared with the Fermi energy of neutrons. 
The baryon density would reach a few times of the nuclear 
density, ~ 10^^ g/cm . Since the neutron density is much 
larger than the proton density, the isospin chemical potential, 
(5/i = (/.i„ — yUp)/2 = (pd — Mu)/2, is finite and large. In rel- 
ativistic mean field (RMF) models, 6fi is calculated to reach 
100 MeV in the neutron star core. Thus we can regard the 
neutron star core matter as asymmetric matter at zero temper- 
ature. 

In Fig.|7] we compare the first order phase transition bound- 
ary in PQM and /3 equilibrium line in RMF at T = B. 
Here we show the boundary for several values of r. RMF 
parameter sets of TMl HI and lOSTY H?! are adopted as 
typical examples. TMl is a model which describes bulk prop- 



erties of normal and neutron rich nuclei as well as the nu- 
clear matter saturation point. lOSTY is an extended version of 
TMl, which includes degrees of freedom of nucleons and hy- 
perons. This comparison shows that for r = 0.2 and 0.3, large 
dfi makes the chiral transition in neutron star crossover In 
lOSTY, hyperons are calculated to appear at /i^ — 1 100 MeV, 
then the transition to quark matter occurs before hyperons 
appear. Since the first order transition generally makes the 
equation of state softer at around the transition density, the 
crossover nature may help to keep the EOS stiff enough and 
to support the heavy neutron stars IBSll . 

We also compare the phase boundary at T = with those 
in flavor SU(3) NJL model results of the neutron star mat- 
ter 1I39II . Since constituent quark mass in Ref. ll39ll is differ- 
ent from that of the present work, we show their results with 
shifted fiB- We find that NJL shows small Sjj. values around 
the transition. This difference mainly comes from the isovec- 
tor coupling with quarks and nucleons. In quark matter, we 
have chosen the vector coupling in the range < r < 0.3. 
In nuclear matter, isovector-vector coupling is chosen to re- 
produce binding energies of neutron rich nuclei, and it corre- 
sponds to 7' ~ 1.0 ^ 1.2. Thus Sfj, is calculated to be larger in 
nuclear matter 



IV. SUMMARY 

We have investigated the QCD phase transition in isospin 
asymmetric matter using the Polyakov loop extended quark 
meson (PQM) model. Specifically, we have discussed isospin 
chemical potential J/^j and quark-vector meson coupling de- 
pendence of the QCD phase boundaries. In PQM, we show 
Sfi reduces the temperature of the QCD critical point, and for 
large Sfi, the critical point is found to disappear We also show 
the finite quark-vector meson coupling shifts the chiral phase 
boundary to higher baryon chemical potential and reduces the 
temperature of the CP. This scenario is in agreement with the 
one obtained within other chiral models ll6l l39ll . 

We have also discussed the order of the chiral phase transi- 
tion in neutron stars from the comparison of the QCD phase 
diagram in PQM and the /3 equilibrium line in RMF. In neu- 
tron stars, 6fi is large, then the temperature of the CP becomes 
lower Therefore the chiral phase transition may be crossover, 
even if the transition in symmetric matter ((5yU = 0) is the first 
order In this study, however, we use (//b, (5yu) values on the 
f3 equilibrium line in neutron star matter calculated with RMF 
models which do not include the QCD phase transition effects. 
In order to discuss the QCD phase transition in compact astro- 
physical phenomena more precisely, we need the EOS which 
includes both baryonic and quark degrees of freedom. 

One may consider the reduction of Tcp shown in this paper 
would contradict to the finite lepton-number chemical poten- 
tial result 1I24I1 . which suggests the insensitivity of Tcp as a 
function of the lepton-number chemical potential. Their re- 
sults correspond to the Si^i range 61^1 < 40 MeV, while we find 
that the shift of Tcp is large in the range Sfi > 50 MeV. Thus 
their results could be consistent with ours. 

The phase diagram structure shown in this article is based 



on the assumption that the s-wave pion condensation is not re- 
alized in dense baryonic matter. This assumption is consistent 
with the functional renormalization group calculation ll23ll and 
s-wave ttN repulsion arguments psl], while the results are 
not in a gree ment with the mean field results of PNJL at fi- 
nite 5^ 0221 ■ As a future work, it is an interesting problem 
to discuss the p-wave pion condensation, the inhomogeneous 
chiral condensate, and the color superconductor phases in the 
three-dimensional thermodynamic variable space, (T, ^, 5fi). 



where Ep = y/p'^ + NP with M = 2a = -4:G{qfqf). Here 
G corresponds to the 4-fermion NJL coupling constant, and 
in agreement with the notation of the main text we have put 
jJ-u ^ 1^ — <^M ^nd fid = i-i + Sii. The last addendum on the 
r.h.s. of the above equation corresponds to the valence quarks 
contributions. The vacuum part is regularized by cutting the 
momentum integral at the scale |p| = A. 

Our strategy is as follows: we perform a Ginzburg-Landau 
expansion of the effective potential. 
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where we have subtracted an irrelevant term which does not 
depend on the condensate. At zero temperature and finite 
chemical potential the coefficients are easily determined from 
an expansion of Eq. (lAlb around ct = 0. We get 



2 4Nc .2 2iVc / 2 2\ 



a^ 
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2-loa 
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Appendix A: Critical point within tlie NJL model at zero 
temperature 

In the main body of this article we have discussed the effect 
of an imbalance of the chemical potentials of u and d quarks 
on the location of the critical point (CP) of the QCD phase 
diagram. Our argument was based mainly on numerical re- 
sults obtained within the PQM model. We found that finite 5fi 
moves CP towards a smaller chemical potential and a lower 
temperature. Therefore, we might expect that a large enough 
Sji causes CP to hit the T = plane, then disappearing from 
the phase diagram. In this Appendix we discuss the same topic 
within the NJL model. We limit ourselves to consider a sys- 
tem of u and d quarks in the chiral limit: this simplifies the 
calculations, and allows to identify unambiguously the loca- 
tion of the chiral phase transition in the phase diagram. Our 
purpose is to show analytically how finite 5fL induces a soft- 
ening of the chiral phase transition at finite /i, pushing the CP 
to lower values of temperature (and baryon chemical poten- 
tial). Eventually, for large enough (5/i the CP hits the T = 
plane. For the purpose of our discussion it is therefore enough 
to consider the system at T = and study the change of the 
order of the chiral phase transition at finite ji. 

The thermodynamic potential of the NJL model at zero tem- 
perature can be written as lUOll 



f} = — - 2N,Nf 

Ct 



(27r)3^ 



dp 

(2^ 



{Ep ~ fif)Q ifif - Ep) ,(A1) 



We notice that ag > causing the potential to be bounded 
from below. As a consequence it is possible to study the phase 
transition studying the signs of the first two coefficients. The 
phase transition is of first (second) order if 04 < (04 > 0). 
At the critical point, where the first and second order transition 
lines meet, one has a2 = 014 = 0. Solving 02 = leads to a 
relationship between /i and 5fi; using the solution of the latter 
in the equation 04 = leads to the critical value of Sfi = Sfic 
at which the CP hits the T = plane, namely 



s^J■c 



2GNcNf 



A2 



1 



(A6) 



Using the standard parameters of the model lllOII we find 
(5/ic ~ 140 MeV. This result shows that finite 5fi changes the 
order of the chiral phase transition at zero temperature and 
finite chemical potential. 

The fact that finite Sfi leads to the softening of the phase 
transition can be grasped from Eq. (IA4t : in fact, for Sfj, -^ fi 
one has 



a4 w 0:4 ((5 /i = 0) 






(A7) 



the above equation shows that (5/Lt ^ makes a^ less negative, 
thus favoring a second order phase transition. The same con- 
clusion can be drawn by using an extended version of the GL 
analysis including derivative terms ll40ll . 
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